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1. Introduction 

The Anderson transition is a continuous zero tempera- 
ture quantum phase transition. Just like continuous ther- 
mal phase transitions, such as the magnetic phase tran- 
sition in spin systems, the critical phenomena of the An- 
derson transition are described by the single parameter 
scaling hypothesis. "'^ It is thought that there are three 
universality classes: orthogonal, unitary, and symplectic. 
In the absence of an applied magnetic field or other effect 
that breaks time reversal symmetry, disordered electron 
systems with spin-orbit coupling belong to the symplec- 
tic universality class. Such systems are the subject of this 
proceeding. 

One of the important quantitative characteristics of 
the Anderson transition is the critical exponent v that 
describes the divergence of the localization length ^ 



^ oc |a: 



(1) 



Here a; is a parameter such as strength of disorder, Fermi 
energy, and Xc is its critical value. For the orthogonal and 
unitary universality classes, v has been determined with 
high precision: v = 1.57 ± 0.02 for the orthogonal class^ 
and v = 1.43 ±0.04 for the unitary class'^ in three dimen- 
sion (3D) with 95% confidence intervals. By comparing 
the numerical estimates of the critical exponent in differ- 
ent models, the universality of these values has also been 
tested.^' ^ 

The value of the critical exponent is expected to be 
helpful in identifying the mechanism that drives the 
metal-insulator transition observed in doped semicon- 
ductors and amorphous materials. The theoretically esti- 
mated value ly can be compared with the experimentally 
measured conductivity exponent ^ 

a{T^O)(x{x~XcT, (2) 

through Wegner's relation^ /i = (c? — 2)i/, d being the 
dimensionality of the system; this predicts ^j, = v in 3D. 
Here a{T ~ 0) is the conductivity at zero temperature. 

In early work, the exponent fi was estimated by ex- 
trapolating the conductivity to zero temperature. More 
recent work avoids this unstable extrapolation by using 
finite temperature scaling.^ This permits not only an es- 
timation of the critical parameters but also a check of 
the validity of the single parameter scaling theory. 

Various estimates of the exponent /i have been ob- 



tained experimentally (see Ref . 6 and references therein) . 
Some claim « 0.5, which seems to violate Chayes et 
aZ.'s inequality^' ^ > 2/3. (Note that according to Weg- 
ner's relation /i = i^.) Stupp et a/.^'^" and Itoh et al.^ 
pointed out that /i « 0.5 is obtained when data that 
may not be sufficiently close to the critical point are an- 
alyzed. When such data are excluded they found that 
/X w 1 is always valid. 

The value /i « 1 is at variance with the theoretical val- 
ues for the 3D orthogonal and unitary classes. To help 
clarify the origin of the discrepancy, a precise theoretical 
estimate for the 3D symplectic class is desirable. Further- 
more, research in spintronics has sparked renewed inter- 
est in the effects of the spin-orbit interaction on quantum 
transport. 

In this proceeding, we report precise analyses of the 
critical phenomena at the Anderson transition in the 
SU(2) model^^'^^ and the three dimensional generaliza- 
tion of the Ando model, ^^'^^ both of which have spin- 
orbit coupling and have symplectic symmetry. We also 
extend our finite size scaling analysis beyond the critical 
region and present a numerical estimate of the (3 func- 
tion. 

2. Model and method 

The Hamiltonian used in this study describes non- 
interacting electrons on a square lattice with nearest 
neighbor hopping 



E 



R{h j)aa'cl^Cja' 



(3) 



We distribute the random potential with box distri- 
bution in the range [—W/2, W/2]. Spin-orbit coupling is 
included in the 2x2 matrix R{i,j). 

In the Ando model, the hopping matrix depends only 
on the direction 



R{i,i-\- e„) 



{n = x,y,z), 



(4) 



where e„ is the unit vector in the n direction and cr„ is 
the Pauli spin matrix. The parameter 9, which is set to 
be = 7r/6 in our simulation, represents the strength of 
spin-orbit coupling. 

In the SU(2) model, the hopping matrix R{i,j) is dis- 
tributed randomly and independently with uniform prob- 
ability on the group SU(2) according to the group invari- 
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ant measure. We parametrize the matrix R{i,j) as 



— e 



cos/3y 
-'-'^^ sin I3i. 



c'"^'^' sin I3tj 



(5) 



and distribute a and 7 with uniform probabihty in the 
range [0, 2tt), and (3 according to the probability density, 
P(/3)d/3 = sin(2/3)d/3 in the range [0,7r/2]. 

We consider a quasi- ID bar of cross section L^. Pe- 
riodic boundary conditions are imposed in the trans- 
verse directions. Wc calculate the quasi-lD localization 
length A with the transfer matrix method^^' and ana- 
lyze the finite size scaling of the dimensionless parameter 
A = X/L. In our simulation we fixed the Fermi energy 
E = and accumulated the data by varying disorder W 
and system size L. 

The single parameter scaling hypothesis implies that 
A obeys 



lnA = F± 



(6) 



Here ^ is the correlation length in the metallic phase and 
the localization length in the localized phase. The sub- 
script ± distinguishes the scaling function in the metallic 
and localized phases. When we analyze the critical phe- 
nomena of the Anderson transition, it is more useful to 
use a different but equivalent form of the single parame- 
ter scaling law 

lnA = Fo(i^/'^^). (7) 

Here V is a smooth function of disorder W that crosses 

zero linearly at the critical point W — Wc- 

Single parameter scaling of A can be described by the 
P functioni4.i8 

We estimate the /? function from the finite size scaling 
analyses of numerical data. The /3 function in the critical 
region is calculated from Fq and the critical exponent v, 
and outside the critical region from F± . 

3. Critical phenomena of the Anderson transi- 
tion 

We calculate A for sizes L = 4,5,6,8,10,12,14 with 
an accuracy of 0.1%. The numerical data and the asso- 
ciated fits are shown in Figs. 1 and 2. In the absence 
of any irrelevant corrections to single parameter scaling 
Eq. (7), plotting A vs W should show the critical disor- 
der as the common crossing point of the data. However, 
the curves for different sizes in Figs. 1 and 2 do not cross 
at a common point, indicating the existence of irrelevant 
corrections to single parameter scaling.^ Such corrections 
to perfect single parameter scaling are usually present in 
a simulation of a finite system. 

To analyse data over the full range of system sizes 
L = [4, 14] , we must take account of corrections to single 
parameter scaling due to an irrelevant variable,^ 



In A = Fo (V'i^/'^) + ^L^Fi [i^L^^"^ 



(9) 



where y is an irrelevant exponent and (pL^ is the corre- 
sponding irrelevant scaling variable. For the purpose of 




Fig. 1. In A as a function of disorder W for the SU(2) model. The 
lines are the fit to the function {nii,ni,mii,mi) = (3,1,2,0). 
Although it may be difficult to confirm the existence of correc- 
tions to single parameter scaling in this figure, finite size scaling 
fit shows that such corrections do exist. 
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Fig. 2. In A as a function of disorder W for the Ando model. The 

lines arc the fit to the function {nii,nj ,m,ji,m,j) = (3,1,2,0). 
There are relatively larger irrelevant corrections to single param- 
eter scaling. 



fitting, we approximate Fq and Fi by their finite order 

expansions 

Fo(x) = In Ac + X + a2x'^ ^ ha„^a;"" (10) 

Fi{x) = 1 + 6ix + --- + 6„,.t"^ (11) 
We also expand ip and </> in terms of the dimensionless 
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Table I. The details of various analyses where irrelevant corrections to single parameter scaling axe considered, and the best fit estimates 
of the critical parameters. is the number of data used in each analysis, Np is the number of parameters, Q is the goodness of fit 

probability. The precision is expressed as 95% confidenec intervals. 



Model 


L 


Nd 




nj 




mj 


Np 


Q 


Wc 


In Ac 


V 


y 


SU(2) 


[4,14] 


279 


3 


1 


2 





10 


0.92 


20.001 ± .017 


-0.616 ± .005 


1.375 ± .016 


-2.5 ± .8 


Ando 


[4,14] 


279 


3 


1 


2 





10 


0.9 


19.099 ± .009 


-0.605 ± .002 


1.360 ± .006 


-3.8 ± .4 



Table II. The details of analyses where corrections to scaling are neglected. The precision is again expressed as 95% confidence intervals. 



Model 


L 


Nd 


riR 


mn 


Np 


Q 


Wc 


In Ac 




V 




SU{2) 


[8,14] 


156 


3 


2 


7 


0.7 


19.984 ± .008 


-0.612 ± 


.002 


1.367 ± 


.007 


Ando 


[10,14] 


115 


3 


2 


7 


0.6 


19.092 ± .013 


-0.603 ± 


.003 


1.361 ± 


.011 



disorder w = {We- W)/Wc, 

tP = ip^w + 1p2w'^ + hV'mHW"" (12) 

= (l>o + (f)iw + ■ ■ ■ + (l)„,w"" . (13) 

The details and results of the analyses are presented in 
Table I. 

By discarding data for smaller systems, an analysis 
of the data that neglects irrelevant corrections becomes 
possible. A reasonable fit is obtained in the SU(2) model 
for L > 8 and in the Ando model for L > 10 (Table II). 
The estimated critical parameters in Table II are in rea- 
sonable agreement with those in Table I. 

The most important point to be drawn from Table I 
and Table II is that the estimates of the exponent for 
the SU(2) model and the Ando model are in almost per- 
fect agreement. As for the irrelevant exponent y, our 
estimates are much less precise and we can not clarify 
whether the values of y for these two models are the 
same or not. 

Figures 1 and 2 show that the movement of the cross- 
ing points in the SU(2) model is smaller than that in the 
Ando model, indicating that the magnitude of irrelevant 
corrections in the SU(2) model is smaller. At the same 
time, we can expect that the spin relaxation length is 
shorter in the SU(2) model because of the random spin- 
orbit coupling. This suggests that the leading irrelevant 
corrections observed for small L in the Ando model might 
be from its relatively larger spin relaxation length. 

4. Estimation of the /3 function 

We estimate the jS function Eq. (8) by accumulating 
more numerical data in the metallic and localized regions 
of the SU(2) model at E = 0. The ranges of disorder and 
system size are W = [4, 19] and L = [8, 14] for the metal- 
lic phase, W = [21,30] and L = [8, 16] for the localized 
phase. The precision of the data for A is 0.3% or 0.6%. 
The data covers the strongly localized to the strongly 
metallic limits. 

We perform five scaling analyses for five regions: 
strongly localized, localized, critical, metallic, and 
strongly metallic. From these analyses we estimate the (3 
function. The methods of analyses are almost the same 
as those used in 2D,-'^^ so wc do not describe them in de- 
tail here. We mention only one difference in the strongly 
metallic region. As discussed in Ref. 18, the asymptotic 




-2-10 1 2 

Fig. 3. Single parameter scaling in the localized region is demon- 
strated. 

value of /3(ln A) may depend on the dimensionality d 

/3{\nA)^d-2 (InA^oo). (14) 

We analyze the data in the strongly metallic region 
supposing that the asymptotic value is one in 3D. For 
large but finite A we speculate that derivations from the 
asymptotic value can be described by an expansion in 
powers of A~^. Stopping at the first order we have 

/3(lnA) = l + | (A»l). (15) 

This corresponds to a linear increase of A with L 

A = a|-6 (A»l). (16) 

Numerical data for large A are well fitted by this form. 

In Figs. 3 and 4 we demonstrate the single parameter 
scaling in the localized and metallic region. We can see 
from these figures that the data for different values of 
disorder W and system size L fall on a common scaling 
curve when expressed as a function of L/^. 

The resulting /3 function is displayed in Fig. 5 as well as 
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Fig. 4. Single parameter scaling in the metallic region is demon- 
strated. 




Fig. 5. The /3 functions for A. The P function in 2D is from 
Ref. 14. 



the P function in 2D, which is reported in Ref. 14. The /3 
functions show non-monotonic behavior. The maximum 
values are Pmax ~ 1-40 in 3D and Pmax ~ 0.15 in 2D. 

5. Summary and Discussion 

In summary, wc analyzed the scaling of the parame- 
ter A in 3D. The critical phenomena of the Anderson 



transition in the SU(2) model and the Ando model are 
analyzed and the critical exponent is estimated with high 
precision. All the estimates of the critical exponent are 
in the range ly = [1.35,1.39]. Our estimate is consis- 
tent with the previous estimates f = 1.3 ± 0.2^^ and 

= 1.36 ± 0.10^° for the symplectic class. We also esti- 
mated the f3 function over the full range from the local- 
ized to the metallic limits. 

Our precise estimate clearly distinguishes the value 
of the critic;al exponent in the symplectic class from 
that for the orthogonal class v = 1.58 ± 0.02.^ On 
the other hand, comparing with the available estimate^ 
v = 1.43 ± 0.04 for the unitary class, the precision is not 
yet sufficient to allow the same statement to be made, al- 
though the theoretical expectation is that the values for 
the unitary and symplectic classes should be different, as 
was clearly demonstrated in 2D.^'^'^'''^^ 

The values of the critical exponent for three univer- 
sality classes of the Anderson transition are inconsistent 
with the exponent /U « 1 measured in experiment. This 
may indicate that the observed transition is not a pure 
Anderson transition and that electron-electron interac- 
tion effects may have to be taken into account to properly 
account for the critical phenomena.^ 
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